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ABSTRACT 

We extend the analysis of Gabor transforms on a Cosmic Microwave Background 



(CMB) temperature map (Hansen, Gorski & Hivon 2002) to polarisation. We study 
the temperature and polarisation power spectra on the cut sky, the so-called pseudo 
power spectra. The transformation kernels relating the full-sky polarisation power 
spectra and the polarisation pseudo power spectra are found to be similar to the 
kernel for the temperature power spectrum. This fact is used to construct a fast power 
spectrum estimation algorithm using the pseudo power spectrum of temperature and 
polarisation as data vectors in a maximum likelihood approach. Using the pseudo 
power spectra as input to the likelihood analysis solves the problem of having to 
invert huge matrices which makes the standard likelihood approach infeasible. 

Key words: methods: data analysis-methods: statistical-techniques: image processing 
cosmology: observations-cosmology: cosmological parameters-polarisation 



1 INTRODUCTION 

Most theories of the early universe predict the temperature and polarisation fluctuations of the CMB to be Gaussian dis- 
tributed. In such models, the angular temperature and polarisation power spectra contain all the information about the 
cosmological parameters which one can determine from observations of the CMB sky. As several combinations of the cosmo- 
logical parameters can give rise to similar temperature power spectra, estimating the polarisation power spectra will break 
the degeneracy and will be of great importance for accurate estimation of cosmological parameters. Also the error bars on 
these parameters can be reduced by exploiting the extra information present in the CMB polarisation power spectra. 



Much effort has been made recently in order to find methods to analyse the CMB temperature 


sower spectrum (Hansen, 


Gorski & Hivon 2002; Oh, Spergel & Hinshaw 199^; Wandeh, Gorski & Hivon 2000; Wandeh 2000; 


Wandelt & Hansen 2001; 


Bond 1995; Bond, Jaffe & Knox 200C; Bartlett et al. 2000; Tegmark 1997; Dore, Knox & Peel 2001; Szapudi et al. 200C; 


Hivon et al. 2002; Balbi et al. 2002). There has been very few publications confronting the even harder task of estimating 



the polarisation power spectra. The framework for analysing the polarisation power spectra has been set in (Zaldarriaga & 
Seljak 1997; Kamionkowski, Kosowsky & Stebbins 1997) but these papers only describe the full likelihood method which is 



far too time consuming also when only considering the temperature power spectrum. In (Tegmark & de Oliveira-Costa 2001) 
a quadratic polarisation power spectrum estimation method was introduced, similar to the one presented in (Pegmark 1997) 
for temperature only. 



In this paper we will extend the method of using the pseudo power spectrum as input to a likelihood estimation procedure 
of the power spectrum as described in ( [Hansen, Gorski fc Hivon 2002 ) (from now on called HGH). We will include the E 
and C mode polarisation pseudo power spectra in the data vector and use techniques similar to those described in HGH to 
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estimate the power spectra. This can be done because the kernels that connect the full sky polarisation power spectra with 
the polarisation pseudo power spectra on an apodised sky are similar to the kernel for the temperature power spectrum. In 
the first part of this paper we will derive the formulae for these kernels and for the polarisation pseudo power spectra and 
discuss their shapes. Then in the second part this will be used for likelihood estimation. 



In this paper the B component polarisation will mostly be neglected. The B polarisation power spectrum is expected to 



be very small and will hardly be detectable by the upcoming MAP or Planck satellite experiments (Jaffe, Kamionkowski & 



Wang 2001). Also the E and B components of polarisation mix on the cut sky as will be discussed in this paper, making the 
B polarisation pseudo spectrum to be dominated by the E component (Lewis, Challinor & Turok 2001; [Chiueh fc Ma 2001 ; 
Tegmark fc de Oliveira-Costa 200% |Bunn et al. 2002| ). 



2 THE GABOR TRANSFORMATION 



In this paper we will use the polarisation power spectra as they are defined in (Zaldarriaga & Seljak 1997). There are three 



polarisation power spectra, Cf, Cf and Cf . These are defined as 
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(4) 

where the coefficients are given in Appendix (^. In the following when we write the full sky power spectrum Ci for 
temperature or polarisation, we will always mean the ensemble averaged power spectrum (Cf). 



As discussed in more detail in HGH, the Gabor transformation of a dataset ( Gabor 1946| ) is just the Fourier transforma- 



tion of the dataset multiplied with a window function called the Gabor window. The window function can be used to cut out 
parts of the dataset in order to study only smaller segments of the dataset at a time. The window can be a top-hat window, 
just setting the unwanted parts of the dataset to zero. Another option is to use a function like a Gaussian to smooth the edges 
between the segment which one wants to study and the parts which are set to zero in order to avoid ringing in the Fourier 
spectrum. In HGH this formalism was extended to the sphere and used for CMB analysis. A disc on the CMB sky was cut 
out, using either a top-hat or a Gaussian window. The Gabor coefficients on this cut sky, called the pseudo power spectrum 
was expressed in terms of the full-sky power spectrum and the kernel connecting the full-sky and the cut-sky power spectrum 
was studied. The aim of this first section is to extend this to polarisation. 

When multiplying the polarisation map with a Gabor window G(n, no) which is axissymmetric about the point no, the 
pseudo spectra on this apodised CMB sky can be written in terms of the full-sky spectra as (Appendix ^ 

{Cf) = Y.Cf,K2{tA') + Y.CfK-2{L,t), (5) 
I' I' 

(Cf) = ^CfK2{l,e!) + ^CfK-2{l,(!), (6) 

I' v 

(Cf) = ^Ci;7^2o(A^'), (7) 
e' 

(8) 

where the kernels 2 and 7^20 are given in Appendix (^. These kernels are expression in terms of the h2{t,l' ,m,m') function 
which is defined similar to the h(l, I' , m, m') function used to express the kernel for the temperature power spectrum in HGH. 

As with the temperature kernels, the polarisation kernels can be evaluated either using the analytical Wigner symbol 



expressions (E5E) and ( E58 ) or faster by direct integration and recursion of the h2{£, (.' ,m,m') functions. The recursion for 
the h2{£, £' ,m,m') functions in Appendix (^) is one of the major results in this paper. This is an extension of the recursion 
for h{ej',m,m') derived in HGH. 
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Figure 1. The kernels K2(i,£') (left plot) and K—2{£,£') (right plot) connecting full and cut sky polarisation power spectra Cf and 
Cf- The left kernel is the one which takes full sky into cut sky and full sky Cf into cut sky . The right kernel is the one 
which mixes the two giving contributions from full sky in cut sky and vice versa. 



In Appendix (|e|) and we show that the polarisation pseudo power spectra are rotationally invariant. For this reason 
we will in the rest of the paper put the centre of the Gabor window on the north pole. This makes the calculations easier 
while keeping the generality of the results. 

In this section we will study the Gabor kernel for a Gaussian and a top-hat Gabor window. The Gaussian is defined as 

G{e) = e-*"/'""") e < Oc, (9) 
G{9) = 9>ec, (10) 

where we will use 5 and 15 degrees FWHM (corresponding to cr = 2.12° and a = 6.38°) and a cut-off angle 9c ~ 3a. We will 
also be comparing with a top-hat window covering the same area on the sky {Oc is the same). The top-hat window is defined 
as 

G{9) = 1 9<9c, (11) 
G{9) = e>9c. (12) 

Studying equation (^) and ^ one sees that the E and B modes are mixing when only a portion of the sky is observed. 
The kernel K2{£,£') is the kernel which takes full sky Cf modes to the pseudo coefficients Cf and similarly for Cf . The 
kernel K-2{l,l') is the one which causes the mixing. In figure (|l]) and (^ we have plotted the kernels K2{£,£') together with 
K-2{£,£') for a 5 and 15 degree FWHM Gaussian Gabor window with dc = 3a. One can see that the diagonal of the K2{£,£') 
kernel is about an order of magnitude larger than the diagonal of the mixing kernel K-2{£,£'). This means that Cf would 
dominate Cf and Cf would dominate Cf provided that the two spectra Cf and Cf were of the same order of magnitude. 
However as discussed before, the Cf are expected to be considerably smaller than Cf in most cosmological models. For the 
E mode this is not a problem as Cf will then hardly be affected. The problem is the B mode which in this case will be 
dominated by the E mode. 



Chiueh & Ma 2001 



The separation of E and B modes of polarisation on t he cut sky was alre ady discussed in ( Lewis, Challinor fc Turok 2001 ; 

[icgmark fc dc Olivcira-Costa 20o"l| ; |Bunn et al. 2002| ). We will in this paper assume that the B mode 



polarisation component is negiible, but note that a possible extension of the power spectrum estimation method outlined here 
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Figure 2. The same as figure (hh for a 15 degree FWHM Gaussian Gabor window. 



to B mode polarisation would be to define 

a+,em = aE,im + iaB.im 
a-,im = as, Im ~ iaB.im, 

with the corresponding pseudo quantities 



J2a-,t',^{H2{e,£' ,m) ^ H.2i£,£' ,m)), 



written in terms of the H2 functions defined in Appendix (^). For the power spectra one gets, 

{a+,fma + ,fm) 



21 + 1 



— {C'£m) + (Cim), 



2£+l 



To get the pseudo power spectra one can use equations ( p^ ) and ( p^ to get 



where the kernels can be written 
1 



K+i£,£') = ^^Y.iH2{£,£,m) + H.2{£,£',m)f, 
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Figure 3. The same as figure for a tophat Gabor window covering the same area on the sky as the Gaussian window in figure 



2£ 



^ J2iH2{e,£, m) - H.2{£, £', m)f. 



(24) 



Separation of E/B modes will not be discussed further in this paper as this is extensively treated in the references above. 

In figure (^) and ^ we have plotted the K2{£,£') and K-2{l,t') kernels for a tophat window covering the same area on 
the sky as the Gaussian windows in figure (|l|) and (^. 

The kernel K2o{t,l') for the cross polarisation power spectrum is shown in figure ^ for a 5 and 15 degree Gaussian 
window and in figure ^ for the corresponding tophat windows. 

As for the temperature kernels, all the polarisation kernels show the same behaviour when changing type and size of 
the window. When going from smaller to larger windows, the diagonals get sharper. Also the tophat kernels have more long 
range correlations than the Gaussian kernels (note that all the plots have the same vertical scale and can be compared directly). 

In figure (^) we have plotted slices of the different kernels at £ = 200 for comparison. The slices are made of the kernels 
for 5 and 15 degree FWHM Gaussian Gabor windows. The first thing to note is that the temperature kernel K {£,£'), the E 
and B kernel K2{£,£') and the temperature-polarisation cross spectrum kernel K2o{£,£') only differ for the far off-diagonal 
elements. At the diagonal their shape and size are the same. For this reason the relation shown in HGH between the width 
of the kernel and the width of the size of the window for the temperature power spectrum is also valid for polarisation. This 
is an important result to be used for the likelihood estimation of the polarisation power spectra in the next section. It shows 
that the number of polarisation pseudo spectrum coefficients to be used in the likelihood analysis should be the same as for 
the likelihood estimation of the temperature power spectrum. 

In figure (^) a similar plot is shown for the corresponding tophat windows. The plot shows that the conclusions made 
for the Gaussian windows are also valid in this case. The shape and size of the three kernels are the same around the 
diagonal. For this reason the results shown for the temperature kernel that the tophat window has larger long range correla- 
tions whereas the Gaussian window has large short range correlations and therefore a wider kernel is also valid for polarisation. 

The kernel K-2{£,£') which mixes the E and B modes on the cut sky is plotted as a dashed line in figures ^ and (^). It 
is much smaller than the three other kernels and the shape seems to differ as well. Note that the height of the mixing kernel 
relative to the other kernels is lower for the 15 degree window than for the 5 degree window. That the size of the mixing 
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Figure 6. Same as figure (0) for the corresponding tophat windows. 



kernel relative to the other kernels is dropping when the size of the window is increasing was to be expected since in the limit 
of full sky coverage the E-B mixing disappears and the mixing kernel must go to zero. 

In figure a slice at ^ = 200 of the temperature kernel and the mixing kernel is shown for the 5 and 15 degree FWHM 
Gaussian Gabor window. The kernels are normalised to one at the peak so that the shapes can be compared. For the Gaussian 
window, the shapes of the kernels still seem to be the same. But the kernels for the corresponding tophat windows shown in 
figure (jl^) does not have a Gaussian shape and differs significantly from the other kernels. 

Since the kernels for the polarisation power spectra have a shape similar to that of the temperature power spectrum the 
effect of a Gabor window on the shape of the power spectrum should also be similar. This can be seen in figure ( pi] ) and 
(p^. The figures show the full sky polarisation power spectra (dashed line) Cf (figure ^l|) and (figure for a standard 
CDM model. In this model the B component of polarisation is zero. On top of the full sky power spectra we have plotted 
the polarisation pseudo power spectra for a 5 and 15 degree Gaussian Gabor window (upper and lower plots respectively) 
normalised so that it can be compared to the full sky spectrum. The pseudo spectra for the corresponding tophat windows 
are plotted as dotted lines. As expected the shape of the polarisation pseudo spectra relative to the full sky spectra is similar 
to that for the temperature spectrum shown in HGH. One difference is that the polarisation pseudo spectra for the Gaussian 
window do not have the characteristic extra peak at low multipole which is seen in the temperature pseudo spectrum. This 
peak in the temperature spectrum arose due to the steep + 1) fall-off of the temperature spectrum at low multipole. 

The polarisation spectra do not have this steep fall-off and for this reason there is no extra peak. 

Because of the mixing of E and B modes there is also a B polarisation component Cf for the pseudo spectrum even 
when the input full sky Cf were zero. This is shown in figure (|l^) where we have plotted the full sky spectrum Cf and the 
pseudo spectra Cf for the 5 and 15 degree FWHM Gaussian Gabor windows and corresponding tophat windows. The pseudo 
spectra are normalised so that they can be compared directly to the full sky spectrum. The dashed lines show the pseudo 
spectra for the Gaussian window. The upper line is for the 5 degree window and the lower line for the 15 degree window. As 
expected the size of the B component is dropping with increasing window size. The Cf for the tophat windows are plotted 
as dotted lines. The shape of the pseudo spectra Cf for the Gaussian windows are roughly following the shape of the full 
sky cf. This could be expected because the mixing kernel K-^2{i,£') for the Gaussian window has a Gaussian shape close to 
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Figure 7. A slice at £ = 200 of tlie icernels combining the full sky and cut sky power spectra. The thick solid line is the kernel K(l,l') 
for the temperature power spectrum, the thin solid line is the kernel K2(J-,f-') for E and B mode polarisation and the dotted line is the 
kernel for the temperature-polarisation cross power spectrum K2o{(-,t.')- AH these kernels go together around the diagonal. They only 
differ for the far off-diagonal elements. The dashed line is the mixing kernel K—2{i,i') which mixes the E and B mode polarisation 
power spectra on the cut sky. This kernel is lower than the other kernels. The upper plot is for a 5 degree Gaussian Gabor window and 
the lower plot for a 15 degree FWHM Gaussian window. 

the diagonal, similar to the other kernels (see figure (§|)). The pseudo spectra Cf for the tophat windows however are much 
smoother due to the much broader kernel (figure ^o|). 

In the same way as for the temperature power spectrum, we have shown that the polarisation pseudo spectra resemble the 
full sky polarisation spectra when the patches on the sky are large enough. This motivates the use of the polarisation pseudo 
power spectra as input to a likelihood estimation of the polarisation power spectra in the same way as for the temperature 
power spectrum showed in HGH. To do likelihood analysis one needs to find theoretical expressions for the correlations 
between diflferent Cf (Z={T,E,C}). 



3 LIKELIHOOD ANALYSIS 

In HGH a Gaussian likelihood ansatz with the pseudo power spectrum as the input data was successfully used to estimate the 
power spectrum. Because of the similarities between the kernels of the polarisation power spectra and the temperature power 
spectrum we can assume that this works for the estimation of the polarisation power spectra as well. We will now show the 
results of some Monte Carlo simulations confirming this assumption. In this section we will assume that the B component of 
polarisation is so small that it can be neglected. We will only concentrate on the T, E and C components as in most standard 
t heories of the early universe the B co mponent will be too small to be measured by the MAP and Planck satellite experiments 
( laffe, Kamionkowski fc Wang 2001 ). 

In figure (|l|) and (|l|) we have plotted the probability distribution of the and from 10000 simulations. The 
probability distribution (histogram) is plotted on top of a Gaussian (dashed line) with mean and FWHM taken from the 
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Figure 8. Same as figure (m) but for tophat windows covering the same area on the sky. 



theoretical expressions derived in Appendix (^. In these simulations we were using a 5 degree FWHM Gaussian Gabor 
window with 6c ~ 3cr. In figure ( p^ and (|l^) we show the results of similar simulations with a 15 degree FWHM Gaussian 
Gabor window. As expected the trend is that the distributions get more and more Gaussian for higher multipoles and for 
bigger windows. For the 15° FWHM window, the distribution is very close to a Gaussian for the multipoles above £ — 50. 

Figure (|l^) and ( p^ show the probability distribution for a tophat window covering the same area on the sky as the 
Gaussian window used in figure (|l^) and (^^. Also this distribution is very close to a Gaussian. 

The previous plots have shown that a Gaussian likelihood ansatz for the polarisation pseudo spectra seems to be a very 
good approximation provided that the window is big enough. As for the temperature spectrum, the approximation is no longer 
valid for the lowest multipoles, but as was shown for the temperature power spectrum, this might only give rise to a very 
small downward bias for the estimates of the lowest multipoles. 



The form of the log-likelihood to minimise is therefore still 
L = d^M"M + lndetM, 



(25) 



where the data vector now consists of the temperature and polarisation power spectra d = {d"^, d^, d"^}. Here the vectors 
are given as 

di = Ci. — {Ce-}, 



where Z = {T,E,B}. Similarly the correlation matrix M will consist of blocks Mzz' defined as 
This structure of the data vector and correlation matrix is shown in figure (po|). 



(26) 
(27) 



For fast likelihood estimation, it is crucial that one can calculate the average pseudo spectra {Cf) and correlation matrix 
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Figure 9. A slice at £ = 200 of the kernel K{£,£') connecting the full sky temperature power spectrum with the cut sky temperature 
power spectrum (solid line) and the kernel K—2{£,£') which is mixing the E and B mode polarisation power spectra on the cut sky. The 
upper plot is for a 5 degree FWHM Gaussian Gabor window and the lower plot for a 15 degree Gaussian window. The kernels are here 
normalised to 1 at the peak ai I = 200 in order to compare the shapes of the kernels. 



M fast. The formalism in HGH which enabled fast calculations of these quantities for the temperature power spectrum is 
extended to polarisation in Appendix (^ (signal) and Appendix (^ (noise). 

In figure (j2l]) we have plotted the signal correlation matrix MJ^F next to the matrix MfgF . A standard CDM power 
spectrum without B mode polarisation was used. The two matrices are very similar. One big difference is that the matrix 
for E mode polarisation is missing the 'wall' at low multipoles present in the temperature matrix. As discussed before this 
is because of the different shapes for the T and E power spectra at low multipoles. The temperature power spectrum drops 
steeply at low I while this is not the case for the E mode polarisation spectrum. 

In figures (|2^) and ( |23| ) the M^F, M^F, MJ^? and M^F matrices are shown. All matrices are diagonally dominant and 
since the values on the diagonals all have the same order of magnitude, they all have to be included in the total matrix M. 



4 RESULTS OF LIKELIHOOD ESTIMATIONS 

The likelihood estimation was carried out in the same way as for the temperature power spectrum. As discussed in HGH, 
when we observe the cut-sky we do not have enough information to estimate the full-sky Ct for all multipoles. One has to 
estimate the Ce in A'''"" bins. Also, the Ci coefficients are not independent on the cut sky, so a limited number A*"'" < A^*'" 
of pseudo spectrum coefficients have to be used as input to the likelihood. How many multipoles depends on the width A£ of 
the Gabor kernel for a given window which is discussed in more detail in HGH. 

The power spectra were estimated in bins defined as 
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Figure 10. This figure is the same as figure M) for tophat windows covering the same area on the sky. 



£{£ + 1) ' 



<l<. 



-6+1 1 



(29) 
(30) 



where £b is the first multipole in bin b. A similar binning does not work for for the temperature-polarisation cross correlation 
power spectrum. The reason for this is the Schwarz inequality Cf' < yJCjCf . During likelihood maximisation one must 
make sure that the estimated value of Cf never exceeds ^JCjCf. The way we solved this problem was to estimate for 
cf / CjCf under the constraint that this value never exceeds 1. So the binning is then 



DfVWD, 



where as before It < 



(31) 
(32) 



512 resolution in Healpix (Gorski, Hivon and Wandelt 1998) and 



As an example we simulated a sky using Nside 
a 10' beam. We added non-uniform noise to the map. A reasonable assumption about the size of the noise deviations for 
polarisation is to take = V^o-f (Zaldarriaga & Seljak 1997). This is what we used in this test, but note that the formalism 



do not require any relation between ctt and ap. The noise level was set so that the signal to noise ration for the temperature 
power spectrum was always well above 1 below the maximum multipole £ = 1024 whereas for the E mode polarisation power 



spectrum it was mostly below 1 (see figure (|2^)). This is close to the values expected for the Planck HFI 143GHz (Bersanelli 
et al. 1996) channel. For the analysis we used a 15 degree Gaussian Gabor window. The result of one single estimation is 
shown in figure (jisl). In this estimation we used A'''" = 100 pseudo spectrum multipoles as input and estimated for A'^*'" — 20 
multipole bins for each of the T, E and C modes. 



To test whether the method is bias or not, we did 60 Monte Carlo simulations. The result of the average of these sim- 
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Figure 11. The windowed polarisation power spectra Cf for a 5 and 15 degree FWHM Gaussian Gabor window cut at 9c = Scr (solid 
line) and for a tophat window covering the same area on the sky (dotted line). All spectra are normalised in such a way that they can be 
compared directly with the full sky spectrum which is shown on each plot as a dashed line. Only in the first plot are all three lines visible. 
In the three last plots, the full sky spectrum and the Gaussian pseudo spectrum (dashed and solid line) are hardly distinguishable. 



ulations is shown in figure (^^. The method seems to be unbiased also for the estimates of the polarisation power spectra. 
Note that the expected noise variance taken from the often used analytic formula for uniform noise (shaded areas on the plot) 
given in (Hivon et al. 2002) and HGH here fails to predict the size of the error bars on the estimates. The expected variance 
taken from the inverse Fisher matrix (dashed lines) fits better with the error bars from Monte Carlo. The reason is that we 
used a noise profile with increasing noise from the centre of the disc and down to the edges, opposite of the Gaussian window. 
This gives the observation with high S/N higher significance in the analysis. This is similar to a result for the temperature 
power spectrum discussed further in HGH. 



5 DISCUSSION 

We have presented a maximum likelihood method to simultaneously estimate the temperature and polarisation power spectra 
from high resolution CMB data in the presence of non-uniform noise and a symmetric Gabor window. An extension of the 
power spectrum estimation method developed in HGH has been made in order to estimate for the polarisation power spectra in 
addition to the temperature power spectrum. In most standard theories for the early universe, the B component polarisation 
will be too small to be observed by the MAP and Planck experiments. For this reason, the method has been tested here under 
the assumption that the B mode polarisation is negligible. In this case the method appears to give unbiased estimates of the 
polarisation power spectra also in the presence of non-uniform noise and a Gabor window. 



The kernels connecting the full sky polarisation power spectra with the cut sky polarisation pseudo power spectra were 
studied and found to be very similar to the kernel for the temperature power spectrum. For this reason the eflect of a cut 
sky and a Gabor window on the polarisation power spectra is similar to the effect on the temperature power spectrum. This 
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Figure 12. Same as figure ( |l]] ) for the temperature-polarisation cross power spectrum ■ 



explains that the method of estimating the power spectrum from the pseudo power spectrum for polarisation was as successful 
as it was for temperature. 

One issue which has not been studied fully here is the inclusion of the B mode polarisation. We demonstrated that the 
E and B mode polarisation power spectra are mixing on the cut sky making detections of the much weaker B component 
difficult. Further work needs to be done in order to include the B component in the likelihood analysis. 

In HGH it was discussed how one can find the noise correlation matrix for temperature using Monte Carlo. This might 
be faster than the analytical approach presented here when the size of the dataset is very huge. With a sufficient number of 
Monte Carlo simulations this was shown to give similar error bars as the analytic treatment. The results for the temperature 
noise matrix is expected to be valid also for the polarisation noise matrices and can be used when the dataset is so big that 
the Monte-Carlo approach is significantly faster, or when correlated noise is present. 

Another extension which was discussed in HGH was the simultaneous analysis of several patches on the CMB sky. This 
extension was shown to work for the temperature power spectrum and we expect that this could work also for polarisation, 
allowing data from several different experiments to be analysed together. In HGH it was shown that with this power spectrum 
estimation method, huge datasets like the ones to be expected from MAP or Planck can be analysed in a reasonable amount 
of time. The most time-consuming step in the method is the construction of the correlation matrix, which in the case of 
polarisation is 9 times longer. This makes also the joint temperature and polarisation power spectrum estimation feasible for 
huge datasets. 
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Figure 13. The full sky power spectrum plotted together with the spectra on the windowed sky. The dashed lines show the 
B spectra for a 5 and 15 degree FWHM Gaussian Gabor window (upper and lower line respectively). The dotted lines are for the 
corresponding tophat windows. The pseudo spectra are normalised so that they can be compared directly with the full sky spectrum. In 
the model used, there was no B polarisation spectrum for the full sky. The shown arise due to the mixing of E and B modes on the 
cut sky only. 
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Figure 14. The probability distribution of Cf" taken from 10000 simulations with a 5° FWHM Gaussian Gabor window truncated at 
9c = 3(T. The variable x is given as x = (Cf — {Cf))/^ {{Cf — {Cf))^). The dashed line is a Gaussian with the theoretical mean and 
standard deviation of the Cf. The plot shows the distribution for i = 50, £ = 200, i = 500, and £ = 800. The probabilities are 
normalised such that the integral over x is 1. 
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APPENDIX A: ROTATION MATRICES 



A spherical function r(n) is rotated by the operator D{aP"f) where a/S'y are the three Euler angles for rotations ( Risbo 199{ ) 
and the inverse rotation is -D(— 7 — P — a). For the spherical harmonic functions, this operator takes the form, 



Y«,„(n') = ^ D^/™(a/37)y«w(n) 



(Al) 



Here is a real coefficient with the foUowing property: 

'^m'm(/5) ^mm' { 0) • 

The D-functions also have the following property: 
(a/37) = 

(a2/3272)-C>^"m(ai/3i7i)i 



(A2) 
(A3) 
(A4) 
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= 50 




= 500 




Figure 15. Same as figure (|l^) for the temperature-polarisation cross spectra 



where (q/37) is the result of the two consecutive rotations (qi/3i7i) and (02/3272) 
The complex conjugate of the rotation matrices can be written as 



D'* , = (-1 



{-m){-mO ■ 



See also Appendix (]^. 



(A5) 



APPENDIX B: SPIN-S HARMONICS 



The spherical harmonic functions Y(m{h) can be generalised to spins harmonics using the rotation matrices in Appendix 
(^. The general definition is 



D-sm{4'2,0, (l>l) 



47r 



sYi,n{0, (l>2)e 



21 + 1 

or in the form which will be mostly used in this paper 



The spin-s harmonics have the orthogonality and completeness relations given by 

dfi sYt,n{n) sYemiii) — Su'Smm' 

sYfrn (n) sVcm (no) = (5(n-no). 



(Bl) 

(B2) 

(B3) 
(B4) 



The complex conjugate of the spin harmonics can be written 
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Figure 16. Same as figure ( |l4| ) for a 15° FWHM Gaussian Gabor window. 



(B5) 



APPENDIX C: SOME WIGNER SYMBOL RELATIONS 

Throughout the paper, the Wigner 3j Symbols wiU be used frequently. Here are some relations for these symbols, which are 
used. The orthogonality relation is, 



E 



I I' I" 

m m' m" j \ m m' M 



= {2l" + ir^5i„L"5ra"M- 



The Wigner 3j Symbols can be represented as an integral of rotation matrices (see Appendix(^)), 



This expression can be reduced to. 



£" \ f £ £' £" 



'"^'"s V mi 7712 ms I \ m[ m'2 m' 



... , (2^ + l)(2^' + l)(2^" + l) £ £' £" \ £ £' £" 

dnYi^{n)Yi,^,{n)Yt„^„{n) = \^ '-^ — '-^ '-[ , // o n n 

" 47r \mmm/\uUU 



(CI) 



(C2) 



(C3) 



APPENDIX D: EXTENSION OF THE RECURRENCE RELATION TO POLARISATION 

For fast calculations of correlation matrices for polarisation, it would be pleasant to have a recurrence relation for 

h2{£,£' ,m,m') = / dnG(n) 2yfm(n) 2y<"m' (n). 



(Dl) 
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Figure 17. Same as figure ( |l^ ) for a 15° FWHM Gaussian Gabor window. 



similar to the one for h(l,£' ,m,m') derived in the Appendix of HGH. Again we simplify the notation by calling the function 
. Separating the spin-2 harmonic one can write 

2Yemi9,(t>)^ 2Afm(cos6i)e~™'*. 

In this way one can write A^,™ as 



(D2) 



dx 2\em{x) 2>^e'm'i^) / d(j)G{9,(l>)e " \ 



(D3) 



where x = cos6'. As before we define 



-'^IttT =y dxx 2>^em{x) 2>^e'm'{x)F^^i{x), (D4) 

where obviously F,nm' (x)* = Fm'm{x) and therefore = {Xp^,)* . The next step is to use a recurrence relation for spin-2 

harmonics 



P(i—l)m 2^(i—l)m{x) , 

where 



Pem 



1 /((^ + l)2-m2)((^+l)2-4) 



£+1 
2rn 



4(^ + 1)2 - 1 



l{l + 1) ' 
In this way one has 

m aH m , Aim 



(D5) 

(D6) 
(D7) 

(D8) 
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Figure 18. Same as figure for a tophat window covering the same area on the sky 



I'm' — Pe'm' A{e' + l)m' 



'^e'm' +P(f'-l)m'^(f'-l),^ 



Subtracting the complex conjugate of equation (09) from equation (D8) the left side is zero and one is left with 



P{l'-l),n' 



.(l'-l)m' , , N .(f'-l)m' 



, <{t'~\)m' -2)m' 

+P(t-l)mAf^f_^^^ - P[e'-2)m'A^^ 



(D9) 

(DIO) 
(Dll) 



This is the final recursion formula. The A^^, elements must be provided before the recurrence is started. Then for each 
(m, m'), set £' = m' + 1 and let £ go from £' and upwards, then set £' — m' + 2 and again let £ go from £' and upwards. 
Continue to the desired size of £' . Note that, in order to get all objects up to ^4^"""°=^ one need to go up to ^ = 2£max for 

each row of £' . This is because of the ^[^^j^ term which demands an object indexed {£ + 1) in the previous £' row. 

To start the recurrence, one can precomputed the A'^"^ factors fast and easily using FFT and a sum over rings on the grid. 
As an example, for the HEALPix grid, we did it the following way. 



Am'm' V ^ \r \ 

^Im — / ^ ^m' m"^lm / ^ 6 



JV,,-1 

\ —2Trij /Nr{m — m ) 



Gr 



(D12) 



j=0 



where the last part is the Fourier transform of the Gabor window, calculated by FFT, r is ring number on the grid and j is 
azimuthal position on each ring. Ring r has A'^,- pixels. 

It turns out that the recurrence can be numerically unstable dependent on the window and multipole, and in order to 
avoid problems we (using double precision numbers) restart the recurrence with a new set of precomputed A^T^, for every 
50th £' row. However for some windows and multipoles the recurrence can run for hundreds of ^-rows without problems. 
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Figure 19. Same as figure ( p^l ) for a tophat window covering the same area on the sky 



APPENDIX E: THE POLARISATION PSEUDO POWER SPECTRUM 

The polarisation spherical harmonic coefficients are defined by means of the tensor spherical harmonics 2Yim{'n) as 

a2,em = J dn2Yil,{n){Q + iU){n), (El) 
a-2,£m = j dn ^2Y;^{n){Q~iU){h), (E2) 
and the inverse transforms are given as 

((3 + i[/)(n) = ^a2,vm' 2Yi,^,{h) (E3) 

t'm' 

iQ-iU){h) = Y,a_2,t'm' -2Yt,rr,'in). (E4) 



It will be advantageous to write these spherical harmonics in terms of the rotation matrices -D^m' defined in Appendix (^. 
Using the formulae in Appendix (pBI) one can write 



2Yt^{n) = ^^DU^{cl>,e,0), 



(E5) 

(E6) 
(E7) 
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Figure 20. The figure shows the structure of the datavector d on the left hand side and the correlation matrix M on the right hand 
side used for joint lil^elihood estimation of temperature and polarisation power spectra. 




Figure 21. The correlation matrices M {£,£') in the figure show the correlations between the temperature pseudo power spectrum 
coefficients and between the E mode polarisation pseudo spectrum coefficients for a 15 degree FWHM Gaussian Gabor window. The left 
plot shows {{CjCj,) - {Cj){Cj,))/{{Cj){Cj,)) and the right plot shows {{CfCf,) - {Cf){Cf,))/{{Cf){Cf,)). A standard CDM power 
spectrum was used to produce the plots. 
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Figure 22. The correlation matrices M{£,£') in the figure show the correlations between the cross-correlation pseudo power spectrum C 
coefficients and between the temperature and cross correlation pseudo spectrum coeflicients C for a 15 degree FWHM Gaussian Gabor 
window. The left plot shows {{CfC^,) - {Cf){Cf,))/i{Cf)(C^,}) and the right plot shows {{CfC^,) - {Cj)(C^,))/i{Cj}{C^,)). A 
standard CDM power spectrum was used to produce the plots. 



The corresponding complex conjugates can be written as (using the relations in Appendix (|A|)) 

2YUh) = ^^l±lD%^ic^,9,0) (E8) 

= y^^(-l)™oL™(0,e,O), (E9) 

_2y;„(n) = ^^^DiU^,9,0) (ElO) 

= y^^(-ir^5-2-™(<^,e,0). (Ell) 

Finally, the power spectrum can be written in terms of a 'divergence free' E component and a 'curl free' B component 

aE,lm ~ —-{ci2,£m+Cl-2,em), (E12) 

aB,im = -i{a2,tm — a-2,im) (E13) 



Now we will define the windowed coefficients aim and Ci for polarisation in an analogous way as for temperature. As in 
HGH we Legendre expand the Gabor Window G{8) which is an axissymmetric function centred at fio, 



G{e) = ^^^-±ig,„P,„(cose) = ^ ff,„y,„,„„(n)y;,„„(no). (E14) 

e"m" 

1 2£" + 1 ^e'' 



= J2 9e"\l-j^DC"{<P,0,0)Yi„m„iho). (E15) 

e"m" 

We define the windowed coefficients aem as 
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NormFJised M({,r) 
(LOG) 





Figure 23. The correlation matrices Mil, I') in the figure show the correlations between the temperature and E mode polarisa- 
tion pseudo spectrum coefficients and between the E mode polarisation and cross correlation pseudo spectrum coefficients C for 
a 15 degree FWHM Gaussian Gabor window. The left plot shows ({CjCf,) - {Cj){Cf,))/({Cj){Cf,)) and the right plot shows 
{{Cfc9) - {Cf){Cf,))/{{Cf){C2)). A standard CDM power spectrum was used to produce the plots. 



a2,em = / dn 2Yf*^{n){Q + iU){h)G{n,no) 



(E16) 



Using the expression for (E3) (Q + iU){h) and writing all 2Yem as D- matrices using expressions (E5), (E6), (E9) and 

^(2^ + l)(2^' + l)(2^" + l) 



(Ell) one gets, 



a2,em ~ 



a-2,e'm' ge" 



(47r)3/2 



-(-ir 



X Yi„^n[ho) i dnD'_2m'{<t>,d,0)D^2-m{4>,0,Q)Dii„{<t),e,Q) 



2^a2,£',n' 2^ gi J—^^ (-1) 



(47r)3/2 ^ ^' 2-K 

I" 



X y^„^/,(no) y dnd7D_2„/(<^,6',7)i32-m(<^,^,7)-C'om"('/',^',7) 
= a2,i'm'h2{i, , m, m, np). 



By using equation (E16), one can also write this as, 
dnG(n, no) 2^/^(11) 2y«'m'(n). 



a-2,en — ^ ] 0-2.1' 
I'm' 



(E17) 

(E18) 

(E19) 

(E20) 
(E21) 

(E22) 



Using the two last expressions, the h2 function can be written in two ways (using relation ( |C2| ) for the last expression), 
dnG(n,no) 2y/,„(n) 2Ye,^,(n) (E23) 



h2{£,£' ,m,m' , no) = 



^ J (2^ + l)(2r + l)(2r + l) ^ 



47r 



(E24) 
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Figure 24. The result of a joint likelihood estimation of the temperature power spectrum (upper plot) and the E (middle plot) and 
C (lower plot) polarisation power spectra. The dotted line shows the full sky average spectrum. The histogram shows the binned input 
pseudo spectrum without noise. The shaded areas around the binned average full sky power spectrum (not shown) show the expected 

deviations from the average using the approximate formula for uniform noise. The bright shaded area shows the cosmic and sample 
variance only whereas the dark shaded area also shows expected variance due to noise. The dots show the estimate with la error bars 
taken from the inverse Fisher matrix. In the analysis a 15 degree FWHM Gaussian Gabor window with a. Oc = 3a cutoff was used. 



X y^/'m" (no) 



-2 2 



£" 




rn 



-rn m" 



(E25) 



As we soon will show, the polarisation pseudo power spectra are rotationally invariant under rotation of the Gabor 
window. For that reason, one can put the centre of the Gabor window on the north pole giving. 



(E26) 
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Figure 25. Same as figure (^) but the dots here are the average of 60 estimates from Monte Carlo simulations. The error bars are the 
average deviations taken from the simulations. The dotted line shows the average full sky spectrum. The shaded areas which are plotted 
around the binned full sky power spectrum (not shown) show the variance taken from the approximate variance formula for uniform 
noise. The dashed lines show the expected variance taken from the inverse Fisher matrix. 



where, 



h2{£,£',m) = h2{£,£',m,m,0) (E27) 
= 2_^gt — (E28) 

X ' (E29) 
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Similarly one gets, 
a-2,im = ^ a^2,e'7nh2{£,£' , -m) 



and 

O.E,ln 



0,3,1 



1 

— ~ — [a2,em + ci~2,e 



aEj'mH2{e,£',m) + i'^aB,i'mH-2{i,i' ,m) 



= ^^aB^i,,„H2{£,e' ,m) - i aE,i'mH-2{i, t' , in) 



(E30) 

(E31) 
(E32) 

(E33) 
(E34) 



Please note that whereas h{l, (! , —m) = h(l, m) a similar relation does not exist for /i2(^, t! , —m). Using the expression 
above, one has that. 



h2ie,e',-m) 



2^91" — (-1) 



I I' I" 
-2 2 



An 

t I' I" 
m — m 



(-1) 



e+i'+e" 



(E35) 



(E36) 



The reason why this is not equal to h2{£,£' ,m) is that the first Wigner symbol is not zero when £ + £' + £" is even, which 
is the case when the whole lower row in the Wigner symbol is 0, as in the case with h{£,£' ,m). It is also obvious from the 
expression (E22). For the scalar case, the relation Y^(_„)(n) — (— l)"'Yf,„(n) ensures that there is no dependency on m in 
h{£,£',m) whereas a similar relation does not exist for the tensor harmonics (but see relation (|Bq)). 



We have further defined, 

H2{£,£',m) = ^{h2{£,£',m)+h2{£,£',-m)) 
H-2{£,£',m) = ^{h2{£,£',m)-h2(£,£',-m)) 
which contrary to h2{£,£' ,m) have an m-symmetry 

H2{£,£',-m) = H2{l,£',m) 
H_2{£,£',-m) = -H^2{t,£',m) 

To find the Ci and (later) the correlation matrices, the following quantities will be needed 



{0'E,lmaE.l' ,m') 



\0'B,lmO'B,l' ,m 



{aE,lmaB,i\m') — 

{0,E,tma*ii ^rn') ~ 
{5-3,(11x5,(1 ^^i) — 



Smm' J2 ^2 ^" ^^)H2{£',£",m) 

- e" 

+ ^Cf„H-2{t,£" ,m)H-2{£' ,£" ,^ 
I" 

5mm' ^C§,H2{£,£",m)H2(£',£" 

- e" 

+ Cf„H-2{t, £", m)H-2{t',£", ' 
I" 

Smm'i Ycf„H2{£,£\m)H-2(£', 
- t" 

+ J2C§'H~2{£,£" ,m)H2{£' ,£" ,m 
I" 

Smm'Ycf„H2{e,£",m)h{e',£",m) 

e" 

-i5mm'Yc'i'"H.2{£,£" ,m)h{£' ,£" ,1 



(E37) 
(E38) 

(E39) 
(E40) 

(E41) 

(E42) 

(E43) 

(E44) 

(E45) 

(E46) 
(E47) 
(E48) 
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For (Cf) one now has, 

(gf) =E ^""27+"'""^ =T.^PMiJ')+Y^C^KMiJ'}. (E49) 

m i' t' 

Where, 

K±2{i, I') = 2^ ^±2(A ^) (E50) 

m 

Using the expression for h2(i,i' ,m), one gets 

K±2{e,e') = i^^^(;ii(^,/,m)±/i2(^,/,m)/i2(^,/,-m)) (E51) 

m 

32;^5 (E52) 

(l±(-l/+*'+^") (E53) 



' /J" 



-220 



327r2 




^ ^' 
-2 2 


L" ' 



)(: 


£' 

— m 





^ \ m -m \ m -m ^ ' 



^ 2 (2f + i)(2r + i) / £ f r X 

= 2.5^ 32;^5 ( _2 2 J 

X (l±(-l)^+*'+*"). (E56) 

Similarly, 

(Cf) = ^Ci?i^2(^,^') + X]Q'^-?^-2(^,^') (E57) 

{C?) = Ycf,K2o{£,e') (E58) 

K2o{lj') = -^^Y^H2ii,i',m)hi£,£',m) (E59) 

m 

^ 2 (2^' + 1) (2^ + 1) ^ ^ £' r \ / £ e' t \ 



2 



167r2 V-2 2 oyVOO 

Before studying those kernels we will first show the rotational invariance. To show that {Cf), {Cf) and {Cf) arc rotationally 
invariant under rotations of the window, one can keeps the dependency on the angle and the expression for the kernel turns 
out to be, 

K±2{l,l' = ^^^Y^(hlie,e',m,m',&o) (E61) 

mm' 

±h2{i,i\m,m\hQ)h2{i^i\ —rn, —m\ho) ) (E62) 



V- V- (2£' + l)V(2£" + l)(2L" + l) 

2^ 2^ 9e"9L" ^ (E63) 



e"L" Tn"M" 

n' a" 



X y,.^.(no)y;.M"(fio)( 2 Q 1( \ ) (E64) 



X 

mm' 



H ( m' -m m" ) ( m' -m M" ) ^^^^^ 



(2£" + l)-l5^„^„5^„jj^„ 

X (l±(-l)^+'^'+'^") (E66) 
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= E^^"^^ -2 2 El^^"-"("'')l (E67) 



47r 

X (1± (-1)*+''+'") (E68) 
= K±2{l,l'). (E69) 
This is independent on the angle no which shows the rotational invariance of the pseudo polarisation power spectra. 



APPENDIX F: ROTATIONAL INVARIANCE 

We now want to show that the pseudo power spectra for polarisation are (as the temperature power spectra) rotationally 

invariant under a common rotation of the sky and window. First, note that the rotation matrices D^^^, are rotating both the 

normal spherical harmonics and the spin-s harmonics. This is easy to show. Assume that one wants to rotate 3y£^(n) with 
the Euler angles (a,/3,7). Using the formula for the normal spherical harmonics one gets, 

.y/°^n) = ^Dl,^(a,/?,7) (Fl) 

= \J ^ Dl,^{a, p, 7)I3i w (0, e, 0) (F2) 

= Di,^{<f>rot,erot,0)^^^ (F3) 

= sYimiArot), (F4) 

where firot is the rotation of the angle n by (a,/3,7). This is clearly general for all spin-s harmonics. Therefore we use the 
method from HGH to show that polarisation pseudo power spectra are rotationally invariant, 

Consider a rotation of the sky and window by the angles (—7 — /3 — a). Then the as,em becomes, 

&lfm = J dfi[6(-7 -P- a)T(n)G(n)] sYeUn). (F5) 
If one makes the inverse rotation of the integration angle n, one can write this as; 

alf^ = / dnr(n)G(n)[£)*(a/37) sY;^{n)], (F6) 



which is just 
-rot 



°L = Y,Dt:'m{af3y) J T(n)G(n) sYL'{&)- (F7) 
The last integral can be identified as the normal as,em- 

a°Li = ^ D^r^'m{0iPl)a.s,im- (F8) 

Thus, 

Y^Dt'^{afi-iraB,im. (F9) 

m' m' 

For CE,im (and analogously for CB,im and Cc,im) one gets 

1 rot rot* /^^1n^ 

— aE,lmaE,im (FIO) 

m 

m' Iea (Fll) 

m m' m" 

^ ^ aE,lm'aE,lrr^"^Dl^'m{0ll3l)Dl^n^{af3')) (F12) 



rirot _ 

Ce,i - Yl 



21 



2£ 



CE,t. (F13) 
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APPENDIX G: THE POLARISATION CORRELATION MATRIX 

To find the correlation matrix M for likelihood estimation of the polarisation power spectra one needs the formulae given 



in equations (E41) to (E4S). As shown there the correlations of the pseudo airn coefficients can be written in terms of the 
h{l,l' ,m) function from IfGH and the H2{i,£' ,rn) and H-2{i,£' ,m) functions. These function can be quickly calculated using 
the important recursion formulae deduced in Appendix (^) and in the Appendix of HGIf. The starting points of these recur- 
sions can also be quickly provided using summations and FFT as explained in HGH. We will now show that the correlation 
function M can be expressed in terms of these functions and for this reason can be calculated quickly. 



(Cl) 



The pseudo power spectra can be written as 

i~T ~T*\ 
2£+l ' 

m 



m 



2£ + l ' 
21 + 1 ■ 



(Gl) 
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To find the correlation function between for polarisation one can follow the same steps as for the temperature correlation 
functions in HGH, and get, 
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where the correlation between are given in equations (E41) to (E48) as sums of h{£,£',m) and H±2(£, £' ,m 



APPENDIX H: POLARISATION WITH NOISE 

Analogously to HGH we will now discuss the noise pseudo power spectra and the noise correlation matrix for polarisation. 
Each pixel in the temperature map has a noise temperature rij and for the polarisation maps we assume and nf to have 
the following properties, 
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We also assume that there is no correlation between noise in the different maps T,Q and U. For the full sky one has, 
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which for this type of noise gives Cf^ = Cf^. 



The pseudo a2,im coefficients can now be found using equations (El) and (E2) we define 
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for an axissymmetric Gabor window G having the value Gj in pixel j. The E and B components are then similarly 
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The correlations between these coefficients are 
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where the last line defines H'2{i,l' ,m,m'). The h'2{£,£' ,m,m') is defined similar to the h'{£,£' ,m,m') function in HGH 
h'2i£,£',m,m') = ^G|(af )^ 2Y/^ 2Yl^„ (H21) 



(H22) 



Note the following relation which was used to obtain equation (|H1 
^G|(af)^ .2Y,^^ .2Yl^, = {~ir+"^'h2{£,£',-m,-m'). 



Again, one can see that when the Gabor window AND noise have azimuthal symmetry this reduces simply to, 

h'2{£,£' ,m,m') = h'2{£,£' ,m). (H23) 
In a similar manner the other aim relations can be found 
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{-ir+"''h2{£,£',~m,-m')] 
= iH'_2(l,l' ,m,m'), (H26) 

where the last line again defines //!_2(i', i?', m, m'). The H±2{1,£' ,m,m') functions which are needed to find the noise correla- 
tion matrices can be quickly calculated using the recursion in appendix (|d|). 

Using these relations one can now find the polarisation pseudo spectra 
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One can further use this to find the noise correlation matrices M^^, g^, , defined as, 
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all others combinations are zero. We then find the total correlation matrix Mzz' ,u' consisting of both signal and noise. As 
for temperature, this is not simply the sum of the correlation matrix for signal and noise, one also gets cross terms. The final 
result is 
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